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Implications of LHC searches for Higgs-portal dark matter 
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The search for the a Standard Model Higgs boson at the LHC is reaching a critical stage as 
the possible mass range for the particle has become extremely narrow and some signal at a mass 
of about 125 GeV is starting to emerge. We study the implications of these LHC Higgs searches 
for Higgs-portal models of dark matter in a rather model independent way. Their impact on the 
cosmological relic density and on the direct detection rates are studied in the context of generic 
scalar, vector and fermionic thermal dark matter particles. Assuming a sufficiently small invisible 
Higgs decay branching ratio, we find that current data, in particular from the XENON experiment, 
essentially exclude fermionic dark matter as well as light, i.e. with masses below « 60 GeV, scalar 
and vector dark matter particles. Possible observation of these particles at the planned upgrade of 
the XENON experiment as well as in collider searches is discussed. 



INTRODUCTION 

The ATLAS and CMS collaborations have recently re- 
ported on the search of the Standard Model (SM) Higgs 
boson with 5 fb^^ data l]. Higgs bosons have been ex- 
cluded in a significant mass range and, ignoring the un- 
likely possibility of a very heavy particle, only the very 
narrow window w 115-130 GeV is now left over. 
There is even a slight excess of events in the data which 
could correspond to a SM like Higgs boson with a mass 
of 125 ± 1 GeV. Although the statistics are not sufficient 
for the experiments to claim discovery, one is tempted 
to take this piece of evidence seriously and analyze its 
consequences. 

In this Letter, we study the implications of these LHC 
results for Higgs-portal models of dark matter (DM) . The 
Higgs sector of the SM enjoys a special status since it al- 
lows for a direct coupling to the hidden sector that is 
renormafizable. Hence, determination of the properties 
of the Higgs boson would allow us to gain information 
about the hidden world. The latter is particularly im- 
portant in the context of dark matter since hidden sector 
particles can be stable and couple very weakly to the 
SM sector, thereby offering a viable dark matter candi- 
date 0. In principle, the Higgs boson could decay into 
light DM particles which escape detection Q. However, 
given the fact that the ATLAS and CMS signal is close to 
what one expects for a Standard Model-like Higgs par- 
ticle, there is little room for invisible decays. In what 
follows, we will assume that 10% is the upper bound on 
the invisible Higgs decay branching ratio, although values 
up to 20% will not significantly change our conclusions. 
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We adopt a model independent approach and study 
generic scenarios in which the Higgs-portal DM is a 
scalar, a vector or a Majorana fermion. We first dis- 
cuss the available constraints on the thermal DM from 
WMAP and current direct detection experiments, and 
show that the fermionic DM case is excluded while in the 
scalar and vector cases, one needs DM particles that are 
heavier than about 60 GeV. We then derive the direct 
DM detection rates to be probed by the XENONIOO- 
upgrade and XENONIT experiments. Finally, we discuss 
the possibility of observing directly or indirectly these 
DM particles in collider experiments and, in particular, 
we determine the rate for the pair production of scalar 
particles at the LHC and a high-energy e+e~ collider. 



THE MODELS 

Following the model independent approach of Ref. Q , 
we consider the three possibilities that dark matter con- 
sists of real scalars S, vectors V or Majorana fermions 
X which interact with the SM fields only through the 
Higgs-portal. The stability of the DM particle is ensured 
by a Z2 parity, whose origin is model-dependent. For 
example, in the vector case it stems from a natural par- 
ity symmetry of abelian gauge sectors with minimal field 
content [5|]. The relevant terms in the Lagrangians are 

ACv = ^^mlV^V^ + ^XviV^Vn^ + ^XhvvH^HV^V^', 

ACf=-\mfXX^\^H^HxX- (1) 

Although in the fermionic case above the Higgs-DM cou- 
pling is not renormafizable, we still include it for com- 
pleteness. The self-interaction terms S"^ in the scalar 
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case and the (V^F^)^ term in the vector case are not es- 
sential for our discussion and we wih ignore them. After 
electroweak symmetry breaking, the neutral component 
of the doublet field H is shifted io H° ^ v + h/y/2 with 
V = 174 GeV and the physical masses of the DM particles 
will be given by 

M| = to| + ]^\hssv^ , 
My = my + ^XhYvv"^ , 



(2) 



In what follows, we summarize the most important for- 
mulas relevant to our study. Related ideas and analyses 
can be found in [^-Q and more recent studies of Higgs- 
portal scenarios have appeared in [HI . 

The relic abundance of the DM particles is obtained 
through the s-channel annihilation via the exchange of 
the Higgs boson. For instance, the annihilation cross 
section into light fermions of mass mfcrm is given by 
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Jn = + 3 X ^fn- There exist different estima- 

tions of this factor and in what follows we will use the 
lattice result — 0.326 [14] as well as the MILC results 
[isj which provide the minimal value Jn — 0.260 and the 
maximal value Jn — 0.629. We note that the most re- 
cent lattice evaluation of the strangeness content of the 
nucleon favors Jn values closer to the lower end of the 
above range. In our numerical analysis, we have taken 
into account these lattice results, which appear more re- 
liable than those extracted from the pion-nucleon cross 
section. 

If the DM particles are light enough, Mdm < ^w/i, 
they will appear as invisible decay products of the Higgs 
boson. For the various cases, the Higgs partial decay 
widths into invisible DM particles are given by 
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where Px = -y^l — 4M|. / m\ . We have adapted the pro- 
gram HDECAY [17] which calculates all Higgs decay 
widths and branching ratios to include invisible decays. 



ASTROPHYSICAL CONSEQUENCES 



where Vr is the DM relative velocity. (The cross section 
for Majorana fermion annihilation was computed in [l2j | 
in a similar framework.) We should note that in our 
numerical analysis, we take into account the full set of 
relevant diagrams and channels, and we have adapted 
the program micrOMEGAs to calculate the relic DM 
density. 

The properties of the dark matter particles can be 
studied in direct detection experiments. The DM inter- 
acts elastically with nuclei through the Higgs boson ex- 
change. The resulting nuclear recoil is then interpreted 
in terms of the DM mass and DM~-nucleon cross sec- 
tion. The spin-independent DM-nucleon interaction can 
be expressed as |3] 
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where m^y is the nucleon mass and /^r parameterizes the 
Higgs-nucleon coupling. The latter subsumes contribu- 
tions of the light quarks {Jl) and heavy quarks (///), 



The first aim of our study is to derive constraints on the 
various DM particles from the WMAP satellite Is^ and 
from the current direct detection experiment XENONIOO 
[igl ]. and to make predictions for future upgrades of the 
latter experiment, assuming that the Higgs boson has a 
mass mh = 125 GeV and is approximately SM-like such 
that its invisible decay branching ratio is smaller than 
10%; we have checked that increasing this fraction to 
20% does not change our results significantly. 

In Fig. [TJ we delineate the viable parameter space 
for the Higgs-portal scalar DM particle. The area be- 
tween the two solid (red) curves satisfies the WMAP 
constraint, with the dip corresponding to resonant DM 
annihilation mediated by the Higgs exchange. We dis- 
play three versions of the XENONIOO direct DM detec- 
tion bound corresponding to the three values of Jn dis- 
cussed above. The dash-dotted (brown) curve around 
the Higgs pole region represents BR'"^ = 10% such that 
the area to the left of this line is excluded by our con- 
straint BR™^ < 10%. The prospects for the upgrade 
of XENONIOO (with a projected sensitivity correspond- 
ing to 60,000 kg-d, 5-30 keV and 45% efficiency) and 
XENONIT are shown by the dotted lines. 

We find that light dark matter, Mdm 60 GeV, vio- 
lates the bound on the invisible Higgs decay branching 
ratio and thus is excluded. This applies in particular to 
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FIG. 1. Scalar Higgs-portal parameter space allowed by WMAP 
(between the solid red curves), XENONIOO and BR'"^ = 10% for 
rrih = 125 GeV. Shown also are the prospects for XENON upgrades. 



the case of scalar DM with a mass of 5-10 GeV consid- 
ered, for instance, in Ref. Q. On the other hand, heavier 
dark matter, particularly for Mum ^ 80 GeV, is allowed 
by both BR"" and XENONIOO. We note that almost the 
entire available parameter space will be probed by the 
XENONIOO upgrade. The exception is a small resonant 
region around 62 GeV, where the Higgs-DM coupling is 
extremely small. 

In the case of vector Higgs-portal DM, the results are 
shown in Fig. [5] and are quite similar to the scalar case. 
WMAP requires the Higgs-DM coupling to be almost 
twice as large as that in the scalar case. This is because 
only opposite polarization states can annihilate through 
the Higgs channel, which reduces the annihilation cross 
section by a factor of 3. The resulting direct detection 
rates are therefore somewhat higher in the vector case. 
Note that for DM masses below mh/2, only very small 
values \hvv<O{10-'^) are allowed if BR'""<10%. 

Similarly, the fcrmion Higgs-portal results are shown 
in Fig. 131 We find no parameter regions satisfying the 
constraints, most notably the XENONIOO bound, and 
this scenario is thus ruled out for A?i///A > 10~^. 

This can also be seen from Fig. SI which displays pre- 
dictions for the spin-independent DM-nucleon cross sec- 
tion (Tsi (based on the lattice /at) subject to the WMAP 
and BR'"" < 10% bounds. The upper band corresponds 
to the fermion Higgs-portal DM and is excluded by 
XENONIOO. On the other hand, scalar and vector DM 
are both allowed for a wide range of masses. Apart from 
a very small region around ^m/j, this parameter space 
win be probed by XENONlOO-upgrade and XENONIT. 
The typical value for the scalar crgj is a few times 10~^ 
pb, whereas tisi for vectors is larger by a factor of 3 which 
accounts for the number of degrees of freedom. 



\ 




1 


Min 

Lattice 

Max 


1 L 
A r.i\ yjiy 1 uu 










XENONUP 


y 








- ■ XENONIT 




WMAP 

1 




— Br= 10% 



50 



100 



150 200 

HM(GeV) 

FIG. 2. Same as Fig. 1 for vector DM particles. 



A,hff/A 



Min 
Lattice 




50 100 150 200 

HM(GeV) 

FIG. 3. Same as in Fig.l for fermion DM; A^/^/A is in GeV""^. 
DARK MATTER PRODUCTION AT COLLIDERS 

The next issue to discuss is how to observe directly the 
Higgs-portal DM particles at high energy colliders. There 
are essentially two ways, depending on the Higgs versus 
DM particle masses. If the DM particles are light enough 
for the invisible Higgs decay to occur, M^m ^ ^"^/i, the 
Higgs cross sections times the branching ratios for the 
visible decays will be altered, providing indirect evidence 
for the invisible decay channel. In the case of the LHC, 
a detailed analysis of this issue has been performed in 
Ref. 7] for instance and we have little to add to it. Nev- 
ertheless, if the invisible Higgs branching ratio is smaller 
than « 10%, its observation would be extremely difficult 
in view of the large QCD uncertainties that affect the 
Higgs production cross sections, in particular in the main 
production channel, the gluon fusion mechanism gg ^ h 
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FIG. 4. Spin independent DM— nucleon cross section versus DM 
mass. The upper band (3) corresponds to fermion DM, the middle 
one (2) to vector DM and the lower one (1) to scalar DM. The 
sohd, dashed and dotted lines represent XENONIOO, XENONIOO 
upgrade and XENONIT sensitivities, respectively. 
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FIG. 5. Scalar DM pair production cross sections at the LHC 
with = 14 TeV as a function of their mass for A^gg = 1 in the 
processes pp -> ZSS, WSS and pp -> W*W*+Z*Z* SSqq. 



[201. I'^ fact, the chances of observing indirectly the in- 
visible Higgs decays are much better at a future e+e^ 
collider. Indeed, it has been shovifn that, at y/s « 500 
GeV coUider with 100 fb~^ data, the Higgs production 
cross sections times the visible decay branc hing fractions 
can be determined at the percent level [U l22j | 

The DM particles could be observed directly by study- 
ing associated Higgs production vifith a vector boson and 
Higgs production in vector boson fusion with the Higgs 
particle decaying invisibly. At the LHC, parton level 
analyses have shown that, although extremely difficult, 
this channel can be probed at the 14 TeV upgrade with 
a sufficiently large amount of data [23| if the fraction of 
invisible decays is significant. A more sophisticated AT- 
LAS analysis has shown that only for branching ratios 
above 30% that a signal can be observed at y/s — 14 
TeV and 10 fb~^ data in the mass range ruh = 100-250 
GeV [13]. Again, at a 500 GeV e+e~ collider, invisible 
decays at the level of a few percent can be observed in 
the process e+e~ — >■ hZ by simply analyzing the recoil of 
the leptonically decaying Z boson [21I, . 

If the DM particles are heavy, Mdm ^ ^"i/i, the situ- 
ation becomes much more difhcult and the only possibil- 
ity to observe them would be via their pair production 
in the continuum through the s-channel exchange of the 
Higgs boson. At the LHC, taking the example of the 
scalar DM particle S, three main processes can be used: 

a) double production with Higgs-strahlung from either a 
W or a Z boson, qq ^ V* ^ VSS with V = W 01 Z, 

b) the WW/ZZ fusion processes which lead to two jets 
and missing energy qq — > V*V*qq — > SSqq and c) the 
gluon-gluon fusion mechanism which is mainly mediated 
by loops of the heavy top quark that couples strongly to 
the Higgs boson, gg h* SS. 

The third process, gg — > 55, leads to only invisible 
particles in the final state, unless some additional jets 



from higher order contributions are present and reduce 
the cross section [2^ and we will ignore it here. For 
the two first processes, following Ref. 26] in which dou- 
ble Higgs production in the SM and its minimal super- 
symmetric extension has been analyzed, we have calcu- 
lated the production cross sections. The exact matrix el- 
ements have been used in the qq — )■ ZSS, WSS processes 
while in vector boson fusion, we have used the longitudi- 
nal vector boson approximations and specialized to the 
WlWl + ZlZl — > 55 case which is expected to provide 
larger rates at the highest energy available at the LHC i.e. 
-/s = 14 TeV (the result obtained in this way is expected 
to approximate the exact result within about a factor of 
two for low scalar masses and very high energies); the 
analytical expressions are given in the Appendix. 

As can be seen from Fig. ?? where the cross sections 
are shown as a function of Mdm for Xhss ~ 1, the 
rates at ^/s = 14 TeV are at the level of 10 fb in the 
WW + ZZ 55 process for Mh < 120 GeV and one or- 
der of magnitude smaller for associated production with 
W and Z bosons. Thus, for both processes, even before 
selection cuts are applied to suppress the backgrounds, 
the rates are small for DM masses of order 100 GeV and 
will require extremely high luminosities to be observed. 

Again, the chances of observing DM pair production 
in the continuum might be higher in the cleaner environ- 
ment of e"'"e~ collisions. The two most important produc- 
tion processes in this context, taking again the example of 
a scalar DM particle, are e+e~ — ZSS that dominates at 



relatively low energies and e'^e — >■ Z* Z* 



-SS 



which becomes important at high energies. The rate for 
WW fusion is one order of magnitude larger but it leads 
to a fully invisible signal, e"'"e~ W*W*vv vvSS. 
Following again Ref. [2g|, we have evaluated the cross 
sections for e^e'^ ZSS at -y/s = 500 GeV (the en- 
ergy range relevant for the ILC) and for Z^Zl — t- 55 
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FIG. 6. Scalar DM pair production cross sections at e+e~ colliders 
as a function of the DM mass for \hSS = 1 in the processes e+e~ 
ZSS at = 500 GeV and ZZ ^ SS at = 3 TeV. 



at = 3 TeV (relevant for the CERN CLIC) and the 
results are shown in Fig|6] as a function of the mass Ms 
for Xhss — 1- One observes that the maximal rate that 
one can obtain is about 10 fb near the Higgs pole in ZSS 
production and which drops quickly with increasing Ms- 
The process ZZ — SS becomes dominant for Ms > 100 
GeV, but the rates are extremely low, below w 0.1 fb. 

The situation should be similar in the case of vector 
and fermion DM and we refrain from discussing it here. 



CONCLUSION 

We have analyzed the implications of the recent LHC 
Higgs results for generic Higgs-portal models of scalar, 
vector and fermionic dark matter particles. Requiring 
the branching ratio for invisible Higgs decay to be less 
than 10%, we find that the DM-nucleon cross section 
for electroweak-size DM masses is predicted to be in the 
range 10~^— 10~* pb in almost all of the parameter space. 
Thus, the entire class of Higgs-portal DM models will be 
probed by the XENONlOO-upgrade and XENONIT di- 
rect detection experiments, which will also be able to 
discriminate between the vector and scalar cases. The 
fermion DM is essentially ruled out by the current data, 
most notably by XENONIOO. Furthermore, we find that 
light Higgs-portal DM A/dm ^ 60 GeV is excluded inde- 
pendently of its nature since it predicts a large invisible 
Higgs decay branching ratio, which should be incompati- 
ble with the production of an SM-like Higgs boson at the 
LHG. Finally, it will be difhcult to observe the DM effects 
by studying Higgs physics at the LHC. Such studies can 
be best performed in Higgs decays at the planned e~^e~ 
colliders. However, the DM particles have pair produc- 
tion cross sections that are too low to be observed at the 
LHC and eventually also at future e+e" colliders unless 
very high luminosities are made available. 

Acknowledgements. The authors would like to thank 



E. Bragina, J.B. de Vivie and M. Kado for discussions 
as well as S. Pukhov for his help in solving technical 
problems related to micrOMEGAs. This work was sup- 
ported by the French ANR TAPDMS ANR-09-JCJC- 
0146 and the Spanish MICINNs Consolider-Ingenio 2010 
Programme under grant Muhi- Dark CSD2009-00064. 



Appendix 



The differential cross section for the pair production 
of two scalar particles in association with a Z boson, 
e+e" — >■ ZSS, after the angular dependence is integrated 
out, can be cast into the form {v ~ 17 A GeV): 



dc7(e^ 



ZSS) 



da:idx2 



(A.l) 

where the electron-Z couplings are defined as = — 1 
and Ve = —1 -I- 4sin^ 9w, 2^1,2 = 2i?i,2/v^ are the scaled 
energies of the two scalar particles, X3 = 2 — xi — X2 is 
the scaled energy of the Z boson; the scaled masses are 
denoted by fj,i = M^ j s. In terms of these variables, the 
coefficient Z may be written as 



Z = - 



1 [iz(x\ \- 8fiz) 



4 (1 - X3 + /iz - 



(A.2) 



The differential cross section has to be integrated over 
the allowed range of the xi,X2 variables; the boundary 
condition is 



2(1 - xi - X2 + 2/15 - ^J■z) + X1X2 



< 1 



(A.3) 



For the cross section at hadron colliders, i.e. for the 
process qq — > ZSS one has to divide the amplitude 
squared given above by a factor 3 to take into account 
color sum/averaging, replace e by g (with — 2/^, Vg — 
21^ — Acq sin^ 9w with and Cq for isospin and electric 
charge) and the cm. energy s by the partonic one s; one 
has then to fold the obtained partonic cross section with 
the quark/antiquark luminosities. The extension to the 
qq — WSS case (with Uq = Vq = v^) is straightforward. 

For the vector boson fusion processes, one calculates 
the cross sections for the 2 2 processes VlVl — SS in 
the equivalent longitudinal vector boson approximation 
and then fold with the Vl spectra to obtain the cross sec- 
tion the entire processes e+e" SSU and qq — >■ qqSS; 
see Ref. [1^ for details. Taking into account only the 
dominant longitudinal vector boson contribution, denot- 
ing by /3y.s the V, S velocities in the cm. frame, s^/^ the 
invariant energy of the VV pair, the corresponding cross 
section of the subprocess VlVl SS reads 



G|M^/ 2 ps 



1 + 13'^ 



1 



Ml) 



(A.4) 
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The result obtained after folding with the vector bo- within about a factor of two for low scalar masses and 
son spectra is expected to approximate the exact result very high energies. 
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